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Abstract — The present study deals with the assessment of 
higher order theory of laminated beams under static 
mechanical loads. The theory has been presented for general 
lay-up of the laminate. The displacement field is expressed in 
terms of only three primary displacement variables by 
satisfying exactly the conditions of zero transverse shear stress 
at the top and bottom and its continuity at layer interfaces. 
The governing equations of motion and boundary conditions 
are derived using virtual work. The number of primary 
displacement unknowns is three, which is independent of the 
number of layers and equal in number to the ones used in the 
first order shear deformation theory. Higher order theory thus 
preserves the computational advantage of an equivalent single 
layer theory. The Third order theory and First order shear 
deformation theory are assessed by comparison with the exact 
two-dimensional elasticity solution of the simply-supported 
beam. A theory is good only if it yields accurate results for all 
kinds of loads and for any lay-up of the beam. For this purpose, 
parametric studies for composite laminates and sandwich 
beams are conducted. 

Index Terms — Composite Beams, First Order Shear 
Deformation Theory, Third Order Theory. 

I. Introduction 

Laminated composite and sandwich beams are used ex- 
tensively in various structures. For their efficient design, a 
good understanding of their deformation characteristics, 
stress distribution, natural frequencies and buckling loads 
under various load conditions are needed. A review of three- 
dimensional (3-D) continuum-based approaches, 2-D theo- 
ries for plates and shells and 1-D theories for beams, along 
with their comparative study for plates under static load- 
ing, has been presented by Saravanos and Heyliger [ 1 ] . Ana- 
lytical 2-D solutions for free vibration [2] are available only 
for simply-supported infinite fiat panels and beams. The 3-D 
finite element analysis [3] results in large problem size which 
may become computationally costly for practical dynamics 
and control problems. Hence, accurate 1-D beam models are 
required without too much loss of accuracy compared to 
the 2-D models. Early works employed elastic beam models 
[4—6] with effective forces and moments due to induced strain 
of actuators. A discrete layer theory (DLT) with layerwise 
approximation of displacements was developed for elastic 
laminated beams with induced actuation strain by Robbins 



and Reddy [7] . Classical laminate theory (CLT) [8] , first order 
shear deformation theory (FSDT) [9] and the refined third 
order theory (TOT) [10,1 1] have been applied without elec- 
tromechanical coupling to hybrid beams. The classical lami- 
nated beam theory (CLT), is an extension of the classical 
beam [13], theory. It neglects the effect of transverse shear 
strain. The CLT ignores the effect of transverse shear defor- 
mation and under predicts deflections and over predicts natu- 
ral frequencies and buckling loads. The first order shear de- 
formation theory (FSDT) based on Timoshenko's beam 
theory assumes linear variation of axial displacement u with 
the thickness coordinate z, across the thickness and the trans- 
verse displacement w to be independent of z. Thus the trans- 
verse shear strain in FSDT is constant across the thickness 
and the transverse shear stress is constant layerwise. Shear 
correction factor is therefore needed to rectify the unrealistic 
variation of the shear strain/stress through the thickness 
which actually even for an isotropic beam is approximated 
parabolic. Higher order shear deformation theories (HSDTs) 
involving higher order terms in the Taylor's expansion of the 
displacements in the thickness coordinate z have been de- 
veloped for laminated composite and sandwich beams. The 
force resultants in FSDT have greater physical meaning than 
the ones used in higher order theories [14, 15] and the num- 
ber of equations of motion is only three. Reddy derived a 
third order variationally consistent theory which satisfies 
the conditions of zero shear stress on the top and bottom 
surfaces of the beam. ID models are used for the beams. In 
equivalent single layer (ESL) theories, global through-the- 
thickness approximations are used for the displacements as 
in CLT, FSDT, and HSDTs. 

ID beam theories in which the displacement field is 
assumed to follow a global variation in the thickness direction 
across all the layers, independent of layup and material 
properties. The deflection w is usually assumed to be constant 
across the thickness and the inplane displacements are 
assumed to follow first the order, third order or higher order 
variations across the thickness. In view of the global variation 
of the displacements along the thickness direction, the 
continuity of transverse shear stress at the layer interfaces 
are not satisfied in these theories. The classical laminate 
theory, the first order shear deformation theory, refined third 
order theory and higher order theories belong to this class. 
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Figure 1. Geometry of composite beam 

Shear correction factors are used in the FSDT, while it is not 
required in the higher order theories. Fig. 1 shows the 
geometry of composite beam. 

II. Governing Equations For Beams 

A. First Order Shear Deformation Theory 

The displacements, strains and stresses in FSDT are 
approximated as: 

u = u D (x , t) + z^pU, t) (1) 
iv = w B (_x, £} (2) 

£ x = ■"d.j + z$d.x (3) 
Yzx = + w ox ( 4 ) 

°x = Qll£*= Qli&d.x + ® 

Without considering the nonlinear terms the nth Fourier 
component, 

KU n - P* (6) 

Where, 

P r = [>i P l PzK 

B. Third Order Theory 

In the third order theory the axial displacement is taken 

as: 

uix. z,t}= u fc, 6) - zwqj + (z - (4z 3 /3^ it )}V r i (7) 
The strain-displacement, 

a, - u (x,t) -zw 0fXX + (z - (4z 3 /3h 2 )yp XiX (8) 
The stresses, 



(9) 



Without considering the nonlinear terms the nth Fourier 
component, 
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Where, 



KU" = P* (10) 
U n = [«„ w„ ^j d ]£, 

III. Results 

A. Types of Beams 

Two highly inhomogeneous simply-supported beams (a) 
and (b) are analysed. The stacking order is mentioned from 
the bottom. Beam (a) is composite beams of material 1 
consisting of four plys of equal thickness 0.25h. Beam (a) 
has symmetric lay-up [0°/W /90° /ft*X The three layer 
sandwich beam (b) has graphite-epoxy faces and a soft core 
with thickness 0. lh /0.8h 10. lh. The orientation 9 k =0 for all 
the plys of beam (b). Beam (a), (b), are shows in Fig. 2 and 
Fig. 3 respectively. 
The material properties are: 
Material 1: 

Y L = 181, F 3 = V 3 = 10.3 GPa 
C 1: = Gjl = 7.17, C 23 = 2,87 GPa 



= 0.25. 



= 0.33 



Face: 



Y t = 131.1, Y z = Y 3 = 6.9 MPa 
G u = 3.588, C S1 = 3.088, = 2.3322 MPa 



Core: 



^= 0.2208, K, = 0.2001, = 2760 MPa 

IS .4 

= 0.99,d, , = u,, = 3 - 10" 5 



G n = 16.56, = 545,1, G 2S =455.4 MPa 



*I2 



90 



Figure 2. Beam (a) 
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Figure 3. Beam (b) 

5. Types of Load 

Consider the following types of load for the analysis: 

1. Sinusoidal load: = ^sin^x/a)- 

2. Uniformly distributed load: p£ = j^. Its non-zero 
Fourier components are (p£} M = 4rpt, fmr for n odd. 

C. Static Response 

Static response is compared for two kinds of load applied 
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Table I. Exact results for static sinusoidal load and for FSDT and TOT 
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TOT 
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-l.U J J4 
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-L.jfff 




Exact 


-2 1393 


-1 1397 


-0.BB16 


-6.3471 


-2.5978 


-1.6469 




FSDT 


0.6312 


06312 


:>.s<:i 


1 1410 


1 1410 


1 1410 


a" JO 5a., 0.5V) 


TOT 


7119 


0.6510 


0.6360 


13139 


1.1834 


1 1510 




Exact 


0.B912 


0.7435 


0.7030 


1.65B0 


1 3519 


1.2725 




FSDT 


-0 6312 


4.6312 


-0 6312 


-1.1410 


-1 1410 


-1 1410 


f JO 5a, -05b.) 


TOT 


417119 


-0.6510 


-0.6360 


-1.3139 


-1.1B34 


-1 1510 




Exact 


-0.8946 


-0.7415 


-0.7023 


-I.673B 


-13529 


-1.2725 




Table II 


Exact results for static uniformly distributed load and for FSDT and TOT 






Beam 


a 


B 




s 




10 


20 


5 


10 


20 




FSDT 


-22391 


-1315S 


-1.0849 


-52737 


-2.5S23 


-1.9095 


w(0.5a.0) 


TOT 


-2.5530 


-1.3975 


-1.1056 


-73339 


-3.1042 


-2.0404 




E.xact 


-2.6748 


-1.4343 


-1.1152 


-7 .3930 


-32521 


-2.077S 




FSDT 


0.8502 


0.8502 


0.8502 


1.5369 


1.5369 


1.5369 


0^(0.53,0.511) 


TOT 


1.0097 


0.8901 


0.8602 


1.8815 


1.6230 


1.5584 




Exact 


1.0711 


0.905S 


0.8641 


1.9655 


1.6441 


1.5637 




FSDT 


-0.8502 


-0.8502 


-O.S502 


-1.5369 


-1.5369 


-1.5369 


0^(0.53,-0.511) 


TOT 


-1.0097 


-0.8901 


-O.S602 


-1.SS15 


-1.6230 


-1.5584 




Exact 


-1.0602 


-0.9031 


-0.S634 


-1.9664 


-1.6443 


-1.5637 


Table III. % Error for FSDT and TOT for static sinusoidal load 
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-10 21 
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TOT 


-20.11 


-12.43 


-9.53 


-20.75 


-12.46 


-9.55 




FSDT 


-29.44 


-14. 86 


-10 12 


-31. S3 


-15.66 


-10.33 


vJOSi, -05b.) 


TOT 


-20.42 




-9.44 


-21.50 


-12.53 


-9.55 



on the top surface of the beam a sinusoidal load p sin(7ix/a), 
a uniformly distributed load p Q . For a thick beam with S = 5 
and for a thin beam with S = 100. Very good convergence of 
the central deflection and stresses is obtained by taking 
Fourier terms from n = lton = N= 159, i.e., for 80 odd terms 
in the Fourier series expansion. The deflection and stresses 
are non-dimensionalised as follows with Y = 6.9 and 10.3 
GPa. 



The central deflection ^(0.5a. 0), axial stresses j (0.5a, 

0.5h) and (0.5a, -0.5h) for the beams (a), (b) under sinusoidal 
static load, obtained by the FSDT and TOT and compare 
with the 2D exact solution, are listed in Table I. for S = 5 ( 
thick beam ), S = 10 ( moderately thick beam) and S = 20 ( thin 
beam ). Similar results for uniformly distributed load are 
presented in Table II. Percentage error for FSDT and TOT for 
static sinusoidal load and for static uniformly distributed load 
are listed in Table III and Table IV. Fig. 4 to Fig. 1 1 shows that 
percentage error varies with respect to h/a for diflection and 
stress. 
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Table IV. % Error for FSDT and TOT for static uniformly distributed load 
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Figure 4. % Error for deflection for test beam (a) under static 
sinusoidal load 



Stirrer 



^((XSc,<l.4h) 




Figure 5. % Error for stress for test beam (a) under static sinusoidal 

load 

IV. Effect Of Angles 

A. Symmetric Laminates 

When ply stacking sequence, material, and geometry (i.e. 
ply thicknesses) are symmetric about the midplane of the 
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Figure 6. % Error for deflection for test beam (b) under static 
sinusoidal load 



Ktttw 



irVQ.SO.ftSJi) 




Figure 7. % Error for stress for test beam (b) under static sinusoidal 

load 

laminate, the laminate is called the symmetric laminate. For a 
symmetric laminate, the upper half through the laminate 
thickness is a mirror image of the lower half. The laminates (- 
45/45/45/-45) = ( - 45 /45 } H and (0/90/90/0) = ( 0/9 ) 9 with 
all layers having the same thickness and material, are example 
of a symmetric cross ply laminate. Now, consider the following 
symmetric laminates with angles (0/0j s ., (0/lS] s .. (0 /30]| SJ 
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(?(0.5a.O) 




Figure 8. % Error for deflection for test beam (a) under static 
uniformly distributed load 



*,(0.S<a,0.Sft) 




Figure 9. % Error for stress for test beam (a) under static uniformly 
distributed load 



P5DT 
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Figure 10. % Error for deflection for test beam (b) under static 
uniformly distributed load 

(0/45) s , and (0/ 60) s . Effect of angles on Deflection and 
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Figure 11.% Error for stress for test beam (b) under static 
uniformly distributed load 

Stress for FSDT and TOT, for S = 5 and for S = 10, for static 
sinusoidal load and for static uniformly distributed load are 
listed in Table V and Table VI. Figs. 12 to 19 shows that effect 
of angles on deflection and stress for static sinusoidal load 
and for uniformly distributed load. 

Conclusions 

A theory is good only if it yields accurate results for all 
kinds of loads and for any lay-up of the beam. Now we can 
compare the results obtained by the FSDT and TOT for the 
static sinusoidal load and the static uniformly distributed 
load are presented in Tables I. and II. For the sinusoidal case, 
the percentage error in deflection and stresses for the 
thickness parameter (S = a/h), S = 5 ( thick beam ), S = 10 ( 
moderately thick beam), and S = 20 ( thin beam ), by TOT less 
than the percentage error obtained by FSDT for both beam 
(a), (b). For the static uniformly distributed load case, the 
percentage error in deflection and stresses for S = 5 ( thick 
beam ), S = 1 ( moderately thick beam) and S = 20 ( thin beam 
) by TOT less than the percentage error obtained by FSDT 
for both beam (a), (b). The TOT yields much more accurate 
results than the FSDT. For the symmetric laminates the effect 
of angles on deflections and stresses for FSDT and TOT for 
static sinusoidal load and for static uniformly distributed load 
are presented in Tables V. and VI. For the sinusoidal case, the 
deflections and stresses for S = 5 and S = 10, are decreases as 
the angle increases and for the static uniformly distributed 
load case, the deflections and stresses for S = 5 and S = 10, 
are also decreases as the angle increases and in both the 
load cases TOT gives much more accurate results than the 
FSDT. 
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